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SUMMARY 


A chemical kinetic mechanism is proposed for the combustion 
of silane-hydrogen^oxygen-nitrogen mixtures in the initial tempera- 
ture range from 800.K to 125QK and pressure range from 0.5 to 1.35 
atm. The mechanism yields results which are in agreement with 
published ignition delay times obtained from shock tube experiments. 

i 

Comparisons between the results obtained using the proposed mechan- 
ism and that set forth in ref. 6 reveal that the former predicts 
appreciably shorter ignition delay times, but a flame blowout 
envelope which is shifted so as to decrease the stable flame region. 
Over much of the thermodynamic range examined, the mechanism pre- 
dicts long reaction times. A three-step global mechanism is pro- 
posed which closely models the ignition phase of SiH^ ~ ^2 ~ a ^ r 
combustion; however, the reaction phase is less well reproduced by 
the global model. The necessity for additional experimental data 
to further assess the proposed models is stressed. 
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INTRODUCTION 


Tite pyrophoric nature of mono-silane (SiH^ , also called silane), 
over a wide range of pressures and temperatures has been well estab- 
lished Ce.g. , refs. 1 and 2). Its pyrophoricity has led to recent 
work using silane as a fuel additive to enhance the ignition, com- 
bustion and flame stabilizing properties of hydrogen (ref. 31) and 
hydrocarbons (ref. 4) burning in air. 

In order to better understand its behavior at the conditions 
encountered in practical devices, such as the combustion of a 
SCRAMJET engine (ref. 3), and to be able to predict a priori its 
behavior at conditions which cannot be readily achieved in the 

laboratory, a reliable chemical 1 kinetic mechanism is required. 

{ 

The first attempt at the formulation of such a mechanism is in 


ref. 3. It was shown subsequently (ref. 5), however, that the 
mechanism in ref. 3 yields predicted ignition delay times which 
are about five times longer than those obtained experimentally in 
shock tube tests. As a result, 'an alternative mechanism was de- 

i 

veloped in ref. 6 which predicted ignition delay times which were 
in agreement with the shock tube data. 

The mechanism in ref. 6 postulates that silane pyrolysis 
SiH 4 -* SiH 2 + H 2 


Cl). 


along with its direct oxidation 


SiH 4 + 0 2 + SiH 3 + H0 2 


( 2 ) 


are the initiating reactions. An examination of the magnitudes 
of the rate constants use in the mechanism of ref. 6 indicates that 
reaction (1) , along with the extremely rapid subsequent oxidation 
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of SiH 2 

S±H 2 + 0 2 + HSiO + OH C3I 

are principally responsible for the disappearance of silane at 
SCRAMJET engine operating conditions . 

Herein, an alternative chemical kinetic mechanism is pro- 
posed which closely models methanei oxidation (.Section 21 . Pre- 
dictions of ignition delay and combustion times ('Section 3) arid 
blowout-limit behavior (.Section 4).: are presented using the pro- 
posed mechanism. Comparisons with: the results predicted using 
the mechanism in ref. 6 reveal substantially different behavior 
over most of the conditions examined. A global chemical kinetic 
mechanism is proposed in Section 5 and the recommendations in 
Section 6 may help to resolve the observed differences between 
the proposed mechanism. 
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SYMBOLS 


kf forward rate coefficient (see Table I) 

m mass flow rate (see Figs. 10-16), g/sec 

p pressure, atm 

T temperature, K 

t time, sec 

V reactor volume (see Figs. 10-16), cm^ 

X percent by volume of silane in the fuel 

R universal gas constant (see Table I) 

<f> equivalence ratio, ratio of fuel (silane plus hydrogen )-to-air 

divided by the stoichiometric fuel-to-air ratio 

Subscripts : 
af adiabatic flame 

ig ignition 

R reaction 

0 initial 

RES residence time 
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CHEMICAL KINETIC MECHANISM 


As was pointed out in ref. 3, "silicon is a meniber of the 
Group IV elements appearing directly below carbon on the periodic 
table." As a result, the silane and methane molecules are 
analogous : 


K[ 

H - C - H 

I 

H 


H 

H - Si - H 
■ 

H 


This chemical-structural similarity may imply that a similar 
chemical kinetic mechanism for the oxidation of the two compounds 
exists. The principal initiating reaction in the oxidation of 
methane is hydrogen abstraction (refs. 7 and 8).: 
ch 4 + M = ch 3 + H + M 

-17 

with a forward rate coefficient of k^ = 2 x 10 exp (-8 8,421.5/RT). . 
By analogy, the reaction 

SiH 4 + M = SiH 3 + H + M ' ~(.i) 

17 

is proposed, having a forward rate coefficient k^ = 2 x 10 x 
exp (-59 , 000/RT) . The activation energy was obtained by scaling 
the value for the equivalent methane reaction using the respective 
C-H and Si-H bond energies . 

In methane mechanisms (refs. 7 and 8), CH 4 + R reactions are 
proposed, in which R = 0, H and OH. Similarly, a sequence of 
SiH 4 + R reactions are proposed here, wherein R = 0,H,0H, SiH 3 and 
H0 3 (Table I) . The reaction with SiH 3 is suggested by the work in 
ref. 9 and the rate coefficient for the SiH 4 + SiH 3 reaction in 
Table I is that for the SiH 4 + CH 3 reaction in ref. 9. The 
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SiH 4 + H.0 2 reaction and rate coefficient are proposed in ref. 6. 

The pyrophoric nature of silane suggests the likelihood of 
the direct attack on SiH^ by 0 2 molecules. In ref. 6, this. is 
treated as a hydrogen abstraction: 

SiH 4 + 0 2 = SiH 3 + H0 2 

It is suggested here that this reaction may, in fact, result in 
the formation of water: 

SiH 4 + 0 2 = SiH 2 0 + H 2 0 (ii) 

As written, reaction (ii) may, in fact, be semi-global in nature, 
being analogous to the thermal decomposition reaction 
SiH 4 SiH 2 + H 2 

in ref. 6 and the reactions 

Cff 3 + 0 2 = CHO + H 2 0 

in ref . 7 and 

HCHO + 0 2 = C0 2 + H 2 0 

in ref. 8. That is, each of the above reactions may embody a 
set of rapidly-occuring elementary reactions whose overall effect, 
relative to the rates of the remaining reactions in the mechanism, 
is summarized by the reaction as written. The pre-exponential 
factor for reaction (ii) was estimated from collision theory and 
the activation energy was obtained from a Semenov-type correla- 
tion Cref . 10 ) . 

The remaining reactions in Table I, along with their rate 
coefficients, were assembled from refs. 3 and 6, with one important 
exception. The pre-exponential factor for the reaction 

SiH 3 + 0 2 = SiH 2 0 + OH (iii) 
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14 3 

is specified in ref. 6 as 8.6 x IQ cm /mole-s, a value which 

exceeds the maximum value predicted by collision theory. The value 

14 3 

used in this study, 1.72 x 10 cm /raole-s, is less than the maxi- 
mum collision theory value. 

Comparisons are made in figs. 1 and 2 between the ignition 
delay times determined experimentally in a shock tube (ref. 5). 
and those predicted by the mechanism in Table I . The agreement 
in figs. 1 and 2 is regarded as quite good. As a result, additional 
studies were undertaken using the mechanism in Table I to assess 
the ignition, combustion and flame stabilization behavior of silane- 
air and silane-hydrogen-air mixtures . 


IGNITION DELAY AND COMBUSTION TIME CALCULATIONS 


The computer program described in ref. 11 was run in its 

constant pressure streamtube mode to determine ignition delay and 

combustion times. Herein, the ignition delay time t. is taken 

to be the time required to go from the initial state of the fuel- 

air mixture to the time when the reacting mixture temperature 

reaches Tq + Q.05 AT, where Tq is the initial temperature. 

At = T - - T r , and T_ £ is the adiabatic flame temperature. The 
at u ctj. 

reaction time t R is defined as the difference between the time 
required to reach Tq +0.95 AT and the ignition delay time. The 
sum of these two times (i.e., the time required to go from the 
initial state Tq to Tq + 0.95 AT) is referred to as the overall 
combustion time . 

Ignition delay times as a function of reciprocal initial 
temperature for various silane-hydrogen fuel mixtures are shown 
in figs. 3 and 4 for one atmosphere and equivalence ratios of 
one and one-half, respectively. The equivalence ratio is the 
overall fuel (silane plus hydrogen). -to-air ratio divided by the 
stoichiometric fuel-to-air ratio. In these figures, silane con- 
centrations range from zero percent (pure H^) to 100 percent 
(by volume) . The substantial reduction of t. , even with as little 
as 2% SiH^ is apparent, particularly at the lower initial tempera- 
tures. Increasing the fuel silane concentration from 2% to 20% 
results in a further decrease of t. by about one order-of-magnitude . 
However, further increases in the fuel silane concentration result 
in relatively small additional reductions of t. . 

The results of ref. 6 are also presented in figs. 3 and' 4 . 
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Except for 2% silane at high temperatures, the mechanism pro- 
posed here predicts appreciably lower ignition delay times. For 
example, for 20% SiH^/8Q% H 2 , Tq = 8Q0K, p = 1 atm, 0 = 1.0, the 
predicted t. is nearly one order-of-magnitude lower than that 
predicted by the mechanism of ref. 6. These substantial dif- 
ferences in the results predicted by the two mechanisms at the 
conditions specified in figs. 3 and 5 will permit an assessment 
to be made of the mechanisms when additional data become available. 

Fig. 5 illustrates the effect of increasing equivalence ratio 
on the ignition delay times of pure silane at one atmosphere. 
Although it is clear that for some value of 0>> | , the pyrolysis 
mechanism represented by eqs . (.1)' and (31 becomes important, the 

results shown in fig. 5 were obtained using the mechanism proposed 
here CTable I). Over the initial temperature range examined 
(800 to 1250K) , the t^ is reduced by about a factor of three 
as 0 is increased from 1 to 4 . However, the mechanism predicts 
that little additional reduction of t^ is achieved by increasing 
the equivalence ratio above 4. 


Fig. 6- shows the effect of pressure on t^ at 0 = 1. For 
all quantities of silane in the fuel (between 2% and 100%). , reduc- 
ing the pressure from one atmosphere to one-half atmosphere about 
doubles the ignition delay time in the initial temperature regime 
examined. 

On the basis of similar studies with pure hydrogen, it was 
anticipated that under many thermodynamic conditions, the reaction 


time would be substantially longer than the ignition delay time; 
hence, t R would become rate-controlling. To assess this, t R -was 
ascertained for many of the cases shown in figs. 3 - 6. Fig. 7 
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presents results for p = 1 atm, 0=1. At these conditions, 
little difference was found in t using 2% and 20% silane in the 
fuel. However, appreciably lower values of t R (ranging from one- 
half to one order-of-magnitude). were calculated for pure silane. 
Comparing these results with those for t^ in fig, 7, the follow- 
ing conclusions may be drawn: 

. For 2% silane, t^ exceeds t R at low temperatures, while 
the reverse is the case at higher temperatures . The cross-over 
occurs at about 890K. 

. For 20% silane, t^ is appreciably shorter than t R , rang- 
ing from about a factor of 6 lower at T Q = 800K to about two 
orders-of-magni tude lower at Tg = 1250K. 

. For pure silane, t. is lower than t_ over the initial 
c lg R 

temperature range examined, ranging from about a factor of 2 
lower at Tg = 800K to about an order-of-magnitude lower at 
Tg = 1250K . 

Values for reaction times at 0=0. 5, p=l atm are shown in 

\ 

fig. 8. In this case, relatively little difference is observed in 
the calculated t values in going from 2% silane to 100% silane. 

As in fig. 7, t. exceeds t R only for 2% silane at temperatures 
between 8Q0K and 895K, while t. is significantly lower than t R 
for the 20% and 100% silane fuels over the initial temperature 
range examined . 

Figure 9 shows t values for 0=1, p = 0.5 atm. Here, the 
2% and 20% silane fuels exhibit virtually the same reaction times, 
while t_, for pure silane is substantially lower. Once again, only 
for the 2% silane fuel is an initial temperature range encountered 
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wherein t.. > t D Cin this case, 8Q0.K < T ft < 850.K1. For the 20% 
and 1Q0% silane cases, t. is appreciably less than t^ over the 
range examined. 

These reaction time results clearly indicate the necessity 
for dealing with both ignition delay and reaction times. Over 
most of the conditions examined, t is substantially greater than 
t^ and, thus, is the rate-controlling phenomenon in the approach 
to the adiabatic flame temperature. 

It is frequently useful to have available a correlation of 
the t^ and t R data discussed above. The form of the correlation 
often used was first proposed in ref . 12. For the data presented 
here, these correlations for the 20% SiH^/80% 0 = 1, cases 

are : 

-9 3 14fi 

t i • = 1.775 x 10 p ' D exp C8254/T Q i sec 
valid to + 2%, and 


1.206 x 10 3 exp (-1.711 x 10 3 T Q ) 
t R = 1.554 • SSC 


valid to + 10% 
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PERFECTLY-STIRRED REACTOR CALCULATIONS 


Engine combustion chambers generally provide regions in 
which the combusting fuel-air mixture may be stabilized subse- 
quent to ignition. These regions generally are the wake-like 
regions behind bluff bodies, struts and rearward-facing steps, 
or the recirculation regions in the vicinity of fuel injection 
locations. An assessment of the flame-holding capabilities of 
various fuel-air mixtures may be made by assuming that these 
regions may be modeled as perfectly-stirred reactors (PSR, see, 
for example, ref. 8). The perfectly-stirred reactor computer 
program used for the studies reported herein was that discussed 
in ref. 13. 

A typical result is shown in fig. 10 which is the standard 
blowout limit representation first proposed in ref. 14; i.e., 
equivalence ratio 0 as a function of the ratio of the mass flow 
rate of reactants into the reactor to the reactor volume, rii/V. 
Since m/V ~ t REE where t RES is the minimum residence time 
required to ensure a stable flame (i.e. , the minimum time an 
element of mass must reside in the reactor to prevent flame 
blowoutl, fig. 10 indicates that the mechanism in ref. 6 predicts 
a larger stable flame region (i.e., a decreased likelihood of 
flame blowout over a wider range of conditions) than does the 
mechanism proposed herein. As was the case with the predicted 
t. values discussed in the previous section, this difference in 
predicted results could provide a basis for assessing the two 
mechanisms should experimental PSR results become available. It 
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should be noted in fig. 10 that, in contrast with such fuels as 
hydrogen and the hydrocarbons, both mechanisms predict a maximum 
iii/V value at a 0 of about 2.5, rather than 0 =1. It is not 
possible to assess, at this time, whether this is physically 
reasonable behavior, or is attributable to either an anomaly in 
the thermochemical data used in the computer program, or in the 
method for determining blowout, or that the reaction mechanism is 
incomplete. 

Blowout limit curves at one atmosphere for pure hydrogen, 

20% silane/80.% hydrogen and 100% silane fuels are shown in figs. 11 
and 12 for initial reactant temperatures of 30 OK and 600K, 
respectively. Similar conclusions may be drawn in both instances: 

. Pure silane flames are most stable at equivalence ratios 
exceeding unity. However, for 0 ~ 1.0 (Tq = 300K) and 0 » 0.8 
IT Q = 600K) , pure hydrogen flames provide greater flame stability. 

In general, the 20% SiH^/80% are less stable than the 
neat fuels, except for the region above 0 1.5 (Tq = 600K) , where 

it is about equivalent to hydrogen, and the region below 0 ~ 0.7, 
where it is about equivalent to silane. 

These conclusions suggest that in engine configurations in 
which 20% SiH^ is supplied to promote hydrogen ignition, pure 
hydrogen should be used to ensure flame stabilization subsequent to 
ignition. 

The results of figs. 11 and 12 are presented in terms of mini- 
mum required residence times (rather than m/V) in figs. 13 and 14. 
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Conclusions similar to those discussed in connection with figs. 11 
and 12 may be drawn from figs. 13 and 14, The reduction in the 
stable flame region resulting from a pressure decrease from 1.0 atm 
to 0.5 atm is shown in figs. 15 and 16. As may be seen in fig. 16, 
the minimum required residence time about doubles in going from 
1.0 atm to 0.5 atm. 
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A GLOBAL SILANE - HYDROGEN - AIR MECHANISM 


Considerations of computer run times and computer storage 
limitations frequently require the use of a global or quasi- 
global chemical kinetic mechanism. This is particularly true when 
developing or using codes describing time-dependent, multi- 
dimensional and/or turbulent flows (e.g. , refs. 8, 15 and 16). 

The following three-step global mechanism is proposed for the 
mechanism proposed herein: 


2 SiH 4 + 0 2 = 2 HSiO + 3H 2 (4) 

2 HSiO + 0 2 = 2 Si0 2 + H 2 (5) 

2 H 2 + °2 = 2 H 2 ° ( 6 ) 

The forward reaction rate coefficients recommended are: 

, _ (0.95234 0 + 4.4762 0 -1 - 4.4286) in 15 

K f , 4 5.644 U ‘ 

P 

T 2 * 30 exp (-19350/RT) 
k. f 5 = 5.0 x 10 10 T 2,80 exp (-14000/RT) 



(2.3829 0 + 76.190 0 _1 - 3.5714) 1n 6 


(- 2 ^-) 0,574 T 3,29 exp (— 16500/RT) 

In the above, X represents the percent by volume of silane 
in the initial fuel mixture. The proposed ranges of validity of 
this, global model are: 


15 


0.5 <_ p 1.0 atm 

0.5 <_ 0 <_ 4.0 

2% < x, < 20% SiH. 

— — 4 

80 OK < T q < 1250K 

Typical results using the proposed global model are presented 
in figs. 17 - 21, wherein the results are compared with those 
obtained using the complete mechanism (Table I) . Over the range 
of conditions examined, the global model well represents the 
ignition phase. However, the subsequent approach to the adiabatic 
flame temperature is less well represented. Studies employing 
global reaction models in large-scale computer codes suggest that 
replication of the ignition period is of principal importance in 
obtaining acceptable results with global models. Hence, use of 
the proposed model should result in an adequate representation of 
the combustion chemistry in such large-scale computer programs. 
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CONCLUDING REMARKS 


A chemical kinetic mechanism has heen developed for the 
oxidation of silane which predicts ignition delay times in agrees 
ment with published shock tube data (ref, 5)_. This mechanism 
differs from one already proposed (ref. 6). primarily in the 
assumed initiating reactions, although other differences, discus- 
sed herein, exist as well. When the two mechanisms are compared, 
the one proposed here generally predicts substantially shorter 
ignition delay times for silane-hydrogen-air combustion over 
most of the range of thermodynamic conditions examined. On the 
other hand, the proposed mechanism predicts a reduced flame- 
stabilizing capability for silane when compared with the mechan- 
ism in ref. 6 . 

These results point to the importance of additional experi- 
ments to provide the data necessary to deduce the mechanism 
which most closely models the data. These experiments include 
additional shock tube tests encompassing a wider range of data 
than that in ref. 5, steady-flow reactor experiments similar to 
those in ref. 17 but without the interfering heterogeneous 
reactions, and stirred reactor experiments similar to those in 
ref. 14 to determine flame stability limits. 

The generally lengthy reaction, times predicted by the 
mechanism should be noted. Since the proposed mechanism is quite 
similar to that in ref. 6 beyond the reaction initiation stage, 
the latter mechanism is likely to predict long reaction times 
as well. 
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TABLE I 

OXIDATION MECHANISM* 

A n E 


H 

2.00E+17 

0 . 

59000. 

H 

7. 76E+11 

0 . 

6980 . 

H20 

2.00E+14 

0 . 

19600. 

H2 

2. 88E+13 

0 . 

2650.0 

OH 

4.10E+12 

0 . 

1580.0 

H20 

8.64E+12 

0 . 

95.000 

H2G2 

3.00E+12 

0 . 

5600. 

OH 

1. 72E+14 

0 . 

11430. 

H 

1. 30E+14 

0 . 

2000. 

H2 

5.00E+12 

0 . 

0.0 

H 

5.00E+16 

0 . 

76600. 

H2 

3. 30E+14 

0 . 

10500. 

OH 

1.80E+13 

0 . 

3080. 

H20 

7.50E.12 

0 

170. 

H02 

3.95E+14 

0.0 

29500. 

H202 

1.00E+12 

0 . 

8000. 

SiO 

5 . OOE+i.4 

0 . 

29000. 

H2 

2.00E+14 

0 . 

0.00 

OH 

1.00E+14 

0 . 

0.00 

H20 

1 . 00E+14 

0 . 

0.00 

H02 

1. 20E+14 

0.0 

3795. 

0 

1.00E+13 

0 . 

6500. 

H 

4.00E+12 

0 . 

5700. 

M 

2.50E+15 

0 .. 

4370. 

0 

7.2E18 

- 1 . 

117908. 

H 

5.5E18 

-1 

103298. 

OH 

5.2E21 

-1.5 

118000. 

M 

2.3E15 

0 . 

-800. 

OH 

1.2E17 

0 . 

45500. 

M 

7.1E18 

- 1 . 

0 . 

OH 

5.8E13 

0. 

18000 . 

H 

2.0E13 . 

0 . 

5166. 

0 

2.2E14 

0 . 

16800. . 

H 

7.5E13 

0 . 

11099. 



TABLE I - Continued 


SILANE OXIDATION MECHANISM* 


Reaction 





A 

n 

E 

H2 

4 

02 

= 

OH 

4 

0H 

1.0E13 

0. 

43000 

H 

+ 

H02 

= 

H2 

+ 

02 

2.4E13 

0. 

695. 

H 

+ 

HO 2 

= 

OH 

+ 

OH 

2.4E14 

0. 

1887. 

H20 

4 

0 

= 

H 

4 

H02 

5.8E12 

.5 

57000 

0 

+ 

HO 2 

= 

OH 

4 

02 

5.0E13 

0. 

1000. 

OH 

+ 

HO 2 

= 

02 

4 

H20 

3.0E13 

0. 

0. 

H02 

+ 

H2 

= 

H 

+ 

H202 

7. 3E12 

0. 

18677 

H202 

+ 

H 

= 

OH 

+ 

H20 

3.2E14 

0. 

8950 

HO 2 

4 

OH 

= 

0 

4 

H202 

5. 2E10 

.5 

21062 

HO 2 

+ 

H20 

= 

OH 

+ 

H202 

2.8E13 

0. 

32785 

H02 

4- 

H02 

= 

H202 

+ 

02 

2.0E12 

0. 

0. 


* The coefficients are for the rate coefficient in the forward direction written 
in the form = AT n exp (-E/RT) , where A is in cmVroole-s for bimolecular 
reactions, ana cmG/mole^-s for termolecular reactions. The activation energy E 
is in cal/mole, R = 1.987 cal/mole-K, and T is the temperature in K. M is any 
third-body molecule. (Note that 2.00E+17 = 2.00 x lO 3 -', etc.) 
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■ Data of ref. 5 
o Proposed mechanism 


5.71 x 10“ a exp(11200/T ) 


\ .98 1.02 1.06 1.10 1. 14 1. 18 1.22 1.26 x 10 

1/T 0 , K' 1 


Variation of ignition delay time with reciprocal temperature for 
mixture I. Silane-oxygen equivalence ratio of 1.0; pressure = 1 






Ignition delay time, ms 


1 1. 18x10 exp(123CX)/T ) 

(ref. 5) 


o Data of ref. 5 
■ Propsed mechanism 


.96 1.00 1.04 1.08 1. 12 1.16 1.20 1.24x10 

l/T , K' 1 
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Fig. 2.- Variation of ignition delay time with reciprocal temperature for 
mixture II. Silane-oxygen equivalence ratio of 0.5; pressure = 
at 
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atmosphere pressure. 
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Fig. 11.- Blowout limit correlations for silane-hydrogen-air mixtures at one 
atmosphere; initial temperature is 300K. 
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Fig. 12.- Blowout limit correlations for silane-hydrogen-air mixtures at one 
atmosphere; initial temperature is 600K. 
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Fig. 15.- The effect of pressure on the blowout limit correlation for 20% SiH^ 
- 80% H 9 fuel; initial temperature is 600K. 
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Fig. 16.- The effect of pressure on the residence time correlation for 20 % 
Si H 4 - 80% H 2 fuel; initial temperature is 600K. 
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Temperature histories for the complete and global combustion 
mechanisms at an equivalence ratio of one-half, one atmosphere, 20% 
SiH/. - 80% Ho fuel. 
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